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ABSTRACT 


This Technical Memorandum is the final deliverable item required by 
Task Assignment 418, Tracker Calibration Using SST Data. It presents 
the results of all studies conducted under this task which were not 
documented previously in Technical Memorandum CSC/TM-75/6110. 

The error analysis studies examine the effects of placing the target 
satellite in an orbit nearly coplanar with the relay satellite and of data 
span length on the accuracy with which the satellite states can be recovered. 
An analysis of error models using actual SST data spans is also included. 
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SECTION 1 - INTRODUCTION 


Satellite-to-satellite tracking (SST) provides a type of orbit observation not 
previously used in orbit determination and data analysis. The observation 
modeling differs from that used in single spacecraft, since it includes the 
dynamics of two satellite orbits, together with the relative geometrical con- 
figuration of the satellites and ground station. Analysis of this type of data 
requires information on the effect of deviations in one satellite orbit on the 
solution obtained for the other satellite orbit, and on the influence of other 
error models on the solutions obtained for both satellites. 

The purpose of Task Assignment 418 is to provide support in the analysis 
of the capability to determine satellite trajectories, system biases, and 
other parameters using SST data. This memorandum presents the results 
of all studies conducted under this task which were not documented pre- 
viously in Technical Memorandum CSC/TM- 75/6110. The studies documented 
in this memorandum include error analyses on simulated and real satellite- 
to-satellite tracking (SST) data, and differential correction processing of 
selected spans of the real SST data. Error analysis studies include; 

• Examination of the effects of placing the target satellite in 

' an orbit nearly coplanar with the relay satellite 

• Examination of the effects of data span length on the accuracy 
of recovery of either or both state vectors when SST data is 
supplemented with trilateration tracking data from the relay 
satellite 

• Analysis of errors modeled in the SST formulation using real 
data from the Applications Technology Satellite-6 (ATS-6) and 
the Geodynamics Experimental Ocean Satellite-3 (GEOS-3) 


These error analysis studies are documented in Section 2. Section 3 presents 
the results of differential correction processir^ of selected spans of real SST 
data to obtain ATS-6 and GEOS-3 state vectors and associated uncertainties. 


SECTION 2 - ERROR ANALYSIS STUDIES 


2. 1 NEARLY COPLANAR ORBITS 

The study described in this section was conducted in order to determine the 
accuracy with which the target satellite orbit can be determined using SST 
data in the case where the tai’get and the relay satellite orbits are coplanar, 
or nearly so. This anticipates situations which may arise in the Tracking 
Data Relay Satellite System (TDRSS), where some of the target satellite 
orbits are expected to be equatorial. The objective of the study was to deter- 
mine the minimum relative inclination of the satellite orbits which permits 
recovery of the target satellite state vector to within a given margin of 
error. 

To that end, 6-hour and 24-hour spans of simulated satellite-to-satellite 
data were generated by the Navigation Analysis Program (NAP) for each of 
eight target satellite orbits having different (but small) Inclination angles 
relative to the relay satellite orbital plane. The simulated relay satellite 
was placed over the Galopagos Islands in a geosynchronous orbit with an 
inclination of approximately one degree relative to the equator. The com- 
panion target satellite orbits were approximately circular with heights of 
840 km, periods of approximately 100 minutes, and inclinations relative to 
tlie relay satellite orbit of 0.0, 0.25, 0.50, 1.0, 2.0, 3.0, 4.0, and 5.0 
degrees. Both range and Doppler data from Rosman were simulated, with 
a 1-minute separation between the data points. 

Starting with a priori state vectors unperturbed from the values widi which 
the dataware originally generated, the simulated data were processed by 
NAP for one Iteration in tlie differential correction mode. The normal matrix 
computed during the differential correction process and the orbit files 
generated by the integrator were passed to the NAP covariance analysis pro- 
gram (NAPCOV), which used this information to compute the contributions 


to solve-for imcertainties arising from the uncertainties in the consider 
parameters and to propagate these quantities in time. 

The analyses were performed by considering uncertainties in the relay 
satellite state vector and in several other parameters. The data weights 
used and the uncertainties assumed for the consider parameters are shown 
in Table 2-1. 

The uncertainties computed for the solve-for parameters were rotated into 
radial, along-traclc, and cross-track coordinates (HLC'^and propagated 
through the data span. The largest uncertainties encountered in the 
propagated span were then tabulated. 

Table 2-2 lists the computed uncertainties in the target satellite position 
{in HLC coordinates) arising from the data noise and from the consider 
parameters for the various relative inclination angles for both the G-hour 
and 24-hour data spans. Note that, for the 6-hour data spans, the radial 
and along-track uncertainties are insensitive to the relative inclination 
angle except for zero degrees, where the uncertainties increase by up to 
40%. However, the cross-track uncertainties show a very sensitive 
dependence on the relative inclination angle, becoming very large when 
the target and relay satellite orbits are precisely coplanar. The sensitivity 
is demonstrated in Figure 2-1, where the maximum uncertainties due to 
consider variables are plotted as a function of the relative inclination angle. 

The computed uncertainties for the 24-hour spans show a similar behavior. 
The radial and along-track uncertainties are insensitive to the relative 
inclination angle, but the cross-track uncertainties increase as the relative 
inclination angle decreases, Figui'e 2-1 illustrates that the cross-track 
uncertainties for the 24-hour spans are smaller and not as sensitive to 
relative inclination angle as those of the 6-hour spans when the orbits 
become coplanar. 


Table 2-1. Uncertainties Assumed for the Measurements 
and for the Consider Parameters 




A; Uncertainties in SST Data (Weighting Factors) 

_8 

Range: 6 x 10 sec = 9 meters (one-way) 

_8 

Doppler: 6x10 sec = 0 , 002 Hz 


B: Uncertainties in Consider Parameters 


Consider satellite position: 

Consider satellite velocity: 

Rosman station location: 

Range bias: 

Doppler bias; 

Solar pressure on ATS-6; 

Coo (coefficient of spherically 
symmetric component of 
gravitational field); 

C^o (coefficient of quadrupole 
term in gravitational field) : 


100 meters (each component) 

1 cm/sec (each component) 

10 meters (each component) 

-7 

10 sec = 15 meters (one-way) 

5 X 10 ^ sec = 0. 002 Hz 

10% 

2 X 10“*^ 

9, 2 X lO"® 
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Table 2-2. Target Satellite Uncertainties as a Function of 

Relative Inclination Angle Between Orbital Planes 
of Target and Relay Satellites 



Relative 

Rlax. Uncertainties Due 

Max. Uncertainties Due 

Significant 

Data 

Inclination 

to Noise 

to Consider Parameters 

Consider 

Span 

Anale tdecr^ 

H 

L 

C 

H 

L 

C 

Parameters 

6 -Hour 

0 

0.08 

0,25 

717.6 

51.4 

255,5 

201000 



0.25 

0.06 

0.20 

21.7 

44.4 

179.7 

4440 

Relay Satel- 


0,5 

0.06 

0.20 

10.7 

44,3 

179.3 

2208 

lite State 


1.0 

0.06 

0.20 

5,02 

44.3 

178.7 

1038 

Vector and 


2.0 

0.06 

0.20 

2.67 

44.3 

179.0 

550.6 

Gravitational 


3.0 

0.06 

0.20 

1.77 

44.6 

180.0 

374.5 

Constant 


4.0 

0.06 

0.20 

1.33 

44.7 

180.0 

282.2 



5.0 

0.06 

0.20 

1.06 

44.7 

180.0 

227.2 


24-Hour 

0 

0.018 

0.064 

14.94 

35.2 

636.7 

5136 



0.25 

0.017 

0.064 

8.04 

36.0 

637.2 

1015 

Relay Satel- 


0.5 

0.017 

|o 064 

3.28 

35.9 

636.6 

492.4 

lite State 
Vector and 


1.0 

0.017 

0.064 

1.62 

35.5 

635.5 

282.0 


2.0 

0,017 

0.064 

0.85 

35. 7 

636.0 

134.2 

Gravitational 


3.0 

0,017 

0.064 

0.57 

35.9 

636.1 

93.3 

Constant 


4.0 

0.017 

0.064 

0.42 

35.9 

635.8 

80.6 


5.0 

0.017 

0.064 

0.34 

36.0 

635.6 

73.0 



Notes: 

1) Results are tabulated for both 6-hour and 24-hour data spans. 

2) Units for radial (H), along-track (L), and cross-track (C) 
uncertainties are meters. 
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Consider Uncertainty (Meters) 





In practice, it is highly unlikely that the target and the relay satellite crbits 
will be precisely coplanar. However, the results presented here demonstrate 
the great advantage to be obtained by ensuring that the relative inclination 
angle is greater than one or two degrees. 

2.2 SIMULATED TRILATERATION DATA 

Error analysis runs were made adding simulated trilateration data to the SST 
data used for the differential correction runs described previously in Section 2 
and Tables 2-2(a) and {b) of Reference 1. The data consisted of 12 hours of 
GEOS-3 SST data taken from Rosman via the ATS-6 relay satellite, augmented 
by equal quantities of coherent-mode data and data relayed by ground transponders 
located at Mojave and GSFC. Tlie coherent-mode and ground-transponder data 
were taken while GEOS-3 was behind the Earth. Error analysis runs were 
made to include 6 hours and 12 hours of this data, thus corresponding to the 
earlier differential correction runs. Analysis was performed in three configur- 
ations: 

a) Solve for the GEOS-3 state vector only, consider uncertainties* in 
the ATS-6 state vector and in several other parameters 

b) Solve for the ATS-6 state vector only, consider uncertainties* in 
the GEOS-3 state vector and in several other parameters 

c) Solve for both the GEOS-3 and the ATS-6 state vectors, consider 
uncertainties* in several other parameters 

For runs in configurations (a) and (b) above, the computed uncertainties in 
the solve-for parameters were rotated to radial, along-track, and cross-track 
(HLC) coordinates and propagated forward through the data span. The results 
listed in Table 2-3 are the maximum values attained by the respective uncer- 
tainties as they were propagated. 

*The uncertainties assumed for the consider parameters and the data weights 
used are shown in Table 2-1. 
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Table 2-3 (a). Results of Error Analysis Runs Using Simulated 
Trilateration and SST Data (Solving for one 
satellite state vector only) 


Solve for GEOS-3 (target) only, consider 14 parameters 


Data 

Span 

Max. Uncertainties Due 
To Noise 

Max. Uncertainties Due 
To Consider Parameters 

Significant 

Consider 


H 

L 

C 

H 

L 

C 

Parameters 

6 -Hour 

0.06 

0.37 

0.15 

62.7 

260.0 

33.0 

ATS Vector and 
Grav. Constant 

12- Hour 

0.04 

0.18 

0.05 

— 

61,9 

429.1 

118.9 

ATS Vector and 
Grav. Constant 


Solve for ATS-6 (relay) only, consider 14 parameters 


Data 

Span 

Max. Uncertainties Due 
To Noise 

Max. Uncertai 
To Consider 1 

nties Due 
’arameters 

Significant 

Consider 


H 

L 

C 

H 

L 

C 

Parameters 

6 -Hour 

0.31 

0.48 

0.52 

2069 

8718 

34100 

GEOS Vector and 
Grav. Constant 

12-Hour 

0.20 

0.41 

0.30 

396G 

24590 

40560 

GEOS Vector and 
Grav. Constant 


Notes: 

1) Units for radial (H), along-track (L), and cross-track (C) 
uncertainties are meters. 

2) Data weights and consider parameter uncertainties are shown in 
Table 2-1. 

3) The uncertainties have been propagated through the data span and 
the maximum values attained by each component are tabulated. 
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Table 2-3{b). Eesults of Error Analysis Runs Using 
Simulated Trilateration and SST Data 
(Solving for both satellite state vectors) 

Solve for both GEOS-3 (target) and ATS-6 (relay), consider 8 parameters 




6-HOUR DATA SPAN 

12-HOUR DATA SPAN 

' 

ELEMENT 

NOISE 
(at Epoch) 

CONSIDEil 
(at Epoch) 

significant 

CONSIDER 

PARAMETERS 

NOISE 
(at Epoch) 

CONSIDER 
(at Epoch) 

SIGNIFICANT 

CONSIDER 

PARAMETERS 

GEOS 

X (meters) 

2. 30 

26.04 

Grav. Constant 

1. 39 

43. 80 

Grav, Constant, 
Range Bias 


Y 

2.07 

22.01 

Grav. Constant 

1.31 

35.97 

Grav. Constant 


Z 

1.81 

27.62 

Grav. Constant, 
Rosman Lat. 

0. 95 

25.03 

Grav. Constant, 
Rosman Lat. 


X (cm/sec) 

0.13 

4.25 

Grav, Constant 

0.07 

4.09 

Grav. Constant, 
RH.*!ge Bias 


Y 

0. 22 

2. 37 

Hosman Lat. , 
Rosman Ht. 

0. 12 

2.69 

Grav. Constant, 
Rosman Lat. 


Z 

0.07 

1.68 

Grav. Constant 

0.03 

1.31 

Gray. Constant 

ATS 

X (meters) 

16.88 

125.0 

Grav. Constant 

10.42 

265.0 

Grav. Constant 


Y 

8.56 

50.3 

Grav. Constant, 
Range Bias 

5. 26 

117.8 

Grav. Constant, 
Range Bias 


Z 

7.36 

211.1 

Grav. constant, 
Ros. Position 

3.60 

164.2 

Grav. Constant, 
Ros, Position 


X (cm/sec) 

0.06 

0.19 

Grav, Const., Rng 
Bias, Ros. Lat. 

0,04 

0.72 

Grav, Constant, 
Range Bias 


Y 

0.12 

0.80 

Grav. Constant 

0.08 

1.85 

Grav. Constant 


Z 

0.06 

1. 63 

Grav. Constant, 
Rosman Lat. 

0,03 

1. 34 

Grav. Constant, 
Range Bias 


Notes: 

1) Data weights and consider parameter uncertainties are shown In Table 2-1. 

2) Uncertainties due to noise and the consider parameters tabulated are those at epoch, which occurred 
approximately 40 minutes before the beginning of the data. 



For runs in configuration (c), limitations in the NAPCOV program prevented 
the simultaneous rotation to HLC coordinates and propagation of two different 
state vectors; hence, the uncertainties for these runs are given in XYZ co- 
ordinates at epoch. 

The following trends are evident in Table 2-3. As expected, when only the 
target satellite is being solved for, the uncertainties arising from noise 
decrease as the length of the data span increases, but those arising from the 
consider variables increase vAtli the length of the data span (for extremely 
short spans of the order of a single pass, the uncertainties arising from both 
the noise and the consider parameters become very large, as shown in 
Reference 1). 

When SST data are included in a solution for the ATS-G state vector only, 
using fixed a priori values for the target satellite state vector, large errors 
may be introduced into the solution unless the measurement uncertainties 
used to weight the SST data are very large. This problem arises because 
discrepancies in the SST measurements, brought about by inaccuracies in the 
target satellite state vector or in dynamic parameters affecting the target 
satellite's orbit, will, if processed on an equal footing with direct tracking, 
produce large and unwanted corrections in the relay satellite's orbit. This 
applies particularly to the Doppler measurement, since the Doppler shift due 
to motion of a geosynchronous relay satellite is very small compared with that 
resulting from the motion of the target. This problem is reflected in the error 
analysis results in Table 2-3 (a), in which the SST and trilateration data are 
given equal weights; the large uncertainties in the ATS-6 state vector arising 
from consider parameters are typical of this effect. If both satellite state 
vectors are adjusted, as in Table 2-3(b), the uncertainties in the solve-for 
parameters become tolerable. Wliile the additional SST data, if properly 
weighted to reflect all the uncertainties affecting the target satellite orbit. 
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should, in. theory, improve the solution, the correct weightings are sufficiently 
large and uncertain to suggest that SST data should be excluded when solving for 
the relay satellite only. 

The major consider parameters are listed in the right-hand column of each 
table. The most prevalent of these are the coefficient of the sphericaily sym- 
metric term in die geopotential (listed as "Grav, Constant" in the tables) 
and the state vector components not being solved for. The range bias term also 
becomes significant for the long data span when both state vectors are being 
solved for. The noise contributions to the solve-for uncertainties are all well 
below those of the consider parameters. 
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2.3 REAL GEOS-3 DATA 

An error analysis study of ATS-G/gEOS- 3 SST data recorded on April 27-28, 
1975 , was conducted in conjunction with, the differential correction data pro- 
cessing of these data described in Section 3. The same observational 
data base was used in both sets of runs, so that the tracking schedules in the 
error analyses correspond to those in the differential correction. The data 
consist of: 

a) Three pairs of consecutive passes of SST relay data from 
GEOS-3 with an interval of approximately 10 hours between 
pairs (for convenience, these six passes are numbered ”1" 
through ”6’' in the discussion to follow). Most of these data 
were Doppler type, with a few range type data points included 
in pass 3. 

b) ATS-6 coherent-mode tracking data from Rosman, with a 
small quantity of data relayed through a ground transponder 
at Santiago. 

c) ATS-6 tracking data from Mojave (only the range data were 
used). 

Tlie data separation for (a) and (b) was 10 seconds, and for (c) was 
1 second. 

The six passes of GEOS-3 data were studied in the following nine segments: 
Pass 1 
Pass 3 
Pass 5 
Passes 1-2. 

Passes 3-4 
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Passes 

5-6 

Passes 

1-4 

Passes 

3-6 

Passes 

1-6 


starting with a priori state vectors determined from the differential 
correction runs (described in Section 3), the data were processed 
by NAP for one iteration in the differential correction mode. 

The normal matrix computed during the differential correction process 
and the orbit files generated by the integrator were passed to the NAP 
covariance analysis program (NAPCOV) to determine the contributions to 
solve-for uncertainties arising from the uncertainties in the consider para- 
meters and to propagate these quantities in time. 

Each span of data was used in three modes; 

1. Solving for both the GEOS-3 and the ATS-6 state vectors 
(while considering uncertainties* in a number of dynamic and 
measurement parameters), using only the SST relay data 

2. Solving for both the GEOS-3 and the ATS-6 state vectors 
(while considering fourteen uncertainties ), using both 
the SST relay data and the adjacent ATS-6 tracking data 

3. Solving for the GEOS-3 state vector only, considering the 
ATS-6 state vector and other uncertainties* 



*The uncertainties are shown in Table 2-4. The Mojave uncertainties are 
not used when solving for GEOS-3 only or when solving for both GEOS-3 
and ATS-6 using only the SST relay data. 
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Table 2-4. Uiicertaiatles Assumed for the Measurements 
and for the Consider Parameters 


A; Uncertainties In SST Data (Weighting Factors) 


Range; 

Doppler: 


-7 

5 X 10 seconds = 75 meters (one-way) 
—6 

2 X 10 seconds = 0. 1 Hz (one-way) 


B; Uncertainties in Consider Parameters 


ATS -6 position;* 

ATS -6 velocity:* 

Rosman station location: 

Rosman range bias; 

Rosman Doppler bias; 

Mojave station location; 

Mojave range bias: 

Solar pressure on ATS- 6 : 

Cqq (coefficient of spherically 
symmetric component of 
gravitational field): 

C^o (coefficient of quadrupole 
term in gravitational field) : 

C 30 (3rd zonal coefficient in 
gravitational field) ; 

C^Q (4th zonal coefficient in 
gravitational field): 


200 meters (each component) 

0.465 cm/sec (X and Y components) 
2.97 cm/sec (Z components) 

10 meters (each component) 

-7 

10 sec = 15 meters (one-way) 

-7 

0. 5 X 10 sec 

10 meters (each component) 

-7 

10 sec = 15 meters (one-way) 

10 % 


2 X lO"*^ 
9.2 X 10“® 

1. 125 X 10 
3. 0 X 10 “® 


*The ATS -6 state vector uncertainties are based on the results of previous 
^ error analysis studies of ATS -6 direct-tracking data. 
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Since NAPCOV was unable to propagate uncertainties in more than one state 
vector at a time, runs in modes 1 and 2 provide only the uncertainties at 
epoch, whereas in mode 3 the uncertainties were rotated to HLC coordinates 
and propagated through the data span. Since some of the consider parameters 
(e.g., the Mojave station coordinates) apply only to the ATS-6 tracking data, 
they contribute to the solve-for uncertainties only in mode 2, 

The results of the runs for modes 1, 2, and 3 are shown in Tables 2-5, 2-6, 
and 2-7, respectively. 

Table 2-5 does not contain results for’ single-pass runs, because a single pass 
of SST relay data by itself is insufficient to solve for both the GEOS~3 and the 
ATS-6 state vectors simultaneously. However, over a time span of two passes, 
the ATS-6 relay satellite position varies relative to the GEOS-3 orbital plane 
sufficiently to allow a solution to be obtained with the SST relay data. However, 
the uncertainties shown in Table 2-4 for the pairs of passes 1-2, 3-4, and 
5-6 are large. Combining several adjacent passes of SST relay data , e.g., 
passes 1-4, 3-6, and 1-6, greatly reduces the large uncertainties due to 
noise, but affords no improvement to the uncertainties due to the consider 
parameters. 

Table 2-5 shows the gx’eat advantage of using the adjacent ATS-6 tracking data 
in conjunction with the SST relay data. For the longer runs, e.g. , passes 1-4, 
3-6, and 1-6, the ATS-6 data reduces the uncertainties due to noise by a factor 
of about 10, and those due to consider parameters by a factor of about 20. 

The two- pass runs, e.g., 1-2, 3-4, and 5-6, show a similar reduction of the 
uncertainties due to both the loise and the consider parameters when the 
adjacent ATS-6 tracking data are included. The very large uncertainties in 
the two- pass data on run span 1-2 arise from the lack of Mojave data in this 
span. Thus, the geometrical advantage of having two different station positions 
was lost. The importance of the Mojave data is further demonstrated by the 
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Table 2-5 GEOS-3 (Target) and ATS-6 (Relay) State Vector Uncertainties 
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Table 2-6. GEOS-3 (Target) and ATS-6 (Relay) State '^ 'ector Uncertainties 
at Epoch Using SST Data and ATS-6 Tracking Data 
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insider parameter uncertainties are shown in Table 2-4. 
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To Consider Paraineters 



ATS-6 state vector, 
grav. coefE. 


ATS-6 state vector, 
grav. coefE. Cqq and 


ATS-6 state vector, 
grav. coeEE. 0^^, 
ATS-6 solar press. 


2013 83.56 


25860 294.0 














ATS-6 state vector, 
grav. coeEE. Cqq 


ATS-6 state vector 
grav coefE. Cnn 


ATS-6 state vector, 
grav. coeEE. *^00 


ATS-6 state vector, 
grav. coeEE. Cqo { 


grav. coeEE. Cqq, ! 
ATS-6 state vector 


ATS-6 state vector, 
grav. coeEE C^^ 


ATS-6 state vector, 
grav. coeEE. C 


1) The maximum values attained by the uncertainties when propagated through 
the data span are tabulated. 

2) The consider parameter uncertainties are those given in Table 2-4, ex- 
cluding the Mojave uncertainties. 

3) Units Eor radial (H), along- track (L), and ci’oss -track (C) uncertainties 
are meters. 










































































the results for pass 1, where, lacking a Mojave component, the data were in- 
sufficient to solve for both the GEOS-3 and the ATS-6 state vectors simultan- 
eously. 

The GEOS-3 uncertainties in pass 3 are much larger than those in pass 5 since 
approximately two-thirds of the SST data points in pass 3 were bad and were 
eliminated fay the editing process. The scarcity of the resulting pass 3 data 
is also reflected in the large uncertainties for the two-pass data span on run 
3-4 shown in Table 2-4. Due to the sparseness of the SST data in pass 3, the 
results shown for run 3-4 in Table 2-4 are based on little more than one pass 
of SST data, which is Insufficient to solve for both the GEOS-3 and the ATS-6 
state vectors simultaneously. 

Table 2-6 contains the results obtained from error analysis runs when solving 
for die GEOS-3 elements only. The noise and consider uncertainties computed 
at epoch in each of these runs were rotated into HLC coordinates and propagated 
through the data span; the values tabulated are the maximum values obtained 
by the uncertainties within the span of propagation. 
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SECTION 3 - DIFFERENTIAL CORRECTION PROCESSING 
OF REAL SST DATA 


The SST data of April 27-28 were used to obtain values for the GEOS-3 and 
ATS-6 state vectors . Differential correction runs were made using NAP on 
the nine segments of the data shown on pages 2-11 and 2-12, In die first run, 
all the ATS-6 direct- tracking data in passes 1-6 were used to solve for die 
ATS-6 state vector only. The resulting state vector* and associated uncer- 
tainties were: 


ATS X - 0. 7091651968036083 x lo"^ meters 

a =170.5 meters 

X 

ATS y = -0.4156110240838270 x 10® 

cr = 31.8 

y 

ATS Z - 0. 2553856432432547 x 10® 

0- = 162. 0 
z 

ATS X - 0. 3030356912466290 x 10® cm/sec 

cr. = 0.27 cm/sec 

X 

ATS Y ^ 0. 5183989881342083 x 10® 

0-^ - 1. 20 

ATS Z = 0. 5191434084728750 x 10^ 

a. 1-22 


With the ATS-6 state vector held constant at this value, the nine segments of 
data were processed by NAP in the differentiai correction mode to solve for 
GEOS-3 only. The results of this study are shown in Table 3-1. In order 
to facilitate comparisons between the solutions obtained, these results are 
specified relative to the GEOS-3 state vector obtained from the segment con- 
taining all six passes of GEOS-3 data, i. e. , run 1-6 in Table 3-1. The state 
vector deviations in Table 3-1 are thus simple mathematical differences be- 
tween die several solutions obtained and the following reference state vector* 
obtained in run 1-6: 

GEOS X -= -0. 12665752957 x lo"^ meters 
GEOS Y - -0. 46002468881 X lo"^ 

GEOS Z = 0. 54157448935 x lo"'^ 

*The epoch is 10^ 36“^ UT on April 27. 
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values listed on page 3-1. No consider variables were included 





















GEOS X ^ -0. 34431286721 x 10® (cm/sec) 

GEOS Y = 0. 53952690741 x 10® 

GEOS Z = 0. 37753037804 x 10® 

It should be noted that the GEOS-3 reference state vector is the final solution 
obtained on one of the data spans, not the a priori value used to process the 
data. The ATS-6 reference state vector was obtained by processing only the 
ATS-6 tracking data, and was subsequently used as an a priori value for all 
the SST data processing. The a priori values used for the GEOS-3 state vector 
were adjusted to obtain convergence for a given data span. In nian 3 cases it 
required two or three runs, adjusting both tlie GEOS-3 a priori vector and the 
data edit criteria, in order to obtain a convergence. 

Although the six solutions which do not include pass 6 are reasonably con- 
sistent, the data on pass 6 appear to have a strong perturbing effect on the 
solutions. These are real data and the reason for this perturbation is un- 
loiown. Since the fit to the data on the ATS-only solution does not become 
appreciably worse at tlie end of the data span, an ATS-6 maneuver is not a 
likely explanation. It is more probable that during the final pass a perturba- 
tion was introduced into either the GEOS-3 trajectory or the phase-locked- 
loop mode of the tracking system. 

The comparatively large computed uncertainties In runs containing pass 3 
reflect the fact that two-thirds of the pass 3 data were edited and not included 
in the solutions . 

An attempt was made to do differential correction runs solving for botli die 
GEOS-3 and the ATS-6 state vectors simultaneously. The various multiple- 
pass data segments were used in tliese runs. In only one of tliese differential 
correction runs, the one utilizing the data in passes 1-4, did the differential 
correction process converge to a solution for GEOS-3 and ATS-6. The re- 
sults, shown in Table 3-2, are expressed as deviations about reference state 
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vectors. The GEOS-3 reference state vector is the same as was used in 
Table 3-1. The ATS-6 reference state vector was chosen to be the ATS-only 
solution shown on page 3-1. The large corrections (see Table 3-2) necessary 
to tills ATS-only solution when processing SST data only indicate an inconsis- 
tency between the ATS-6 tracking data and the SST data. 
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APPENDIX A - MATHEMATICAL METHOD USED IN ERROR ANALYSES 


The notation in tliis section is taken from Cliapter 8 of Reference 2. In 
particular, the following are defined: 

= a priori value of vector of solve -for parameters 

= a priori value of vector of consider parameters 

W = diagonal measurement weighting matrix 

f (x, z) = measurement modeling algorithm as a function of solve-for and 
consider parameters 


F 


E 



partial derivative matrix with respect to 
solve-for parameters 


5z 


"o> 


partial derivative matrix with respect to 
consider parameters 



a priori (Bayesian) covariance matrix for solve-for parameters 


P = a prion covariance matrix for consider parameters 
Az 

0 


The method used by NAP and NAPCOV ignores the correlation matrix between 
solve-for and consider parameters, and in ail studies reported in this document 
the Bayesian term P^ was also omitted. With these simplifications, Equa- 
tion ^-39) of Reference 2 reduces to: 

p =i;^rF'^WEP, e'^WF-!- I 
Ax L Az J 

0 

where i/f = (f"^ W F)"^ 

P , = f"^ W E P^ WF + ij/ 

Ax Az 
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The first term of this expression gives the covariance of the solve-for param- 

ettors arising from the consider parameter covariance P ; the second term 

Az 

0 

reflects only the data noise. The square roots of the diagonal elements of 
these matrices provide the uncertainties in solve-for parameters due to con- 
sider parameters and noise, which are listed in the body of this report. Re- 
placing the complete diagonal matrix P by a matrix containing a nonzero 

f\7u 

element for only one consider parameter allows tlie effects of individual 
consider parameters to be estimated; in this way the major contributoi’s to 
the overall uncertainties are identified. 
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